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Abstract 

A  mathematical  model  for  simulating  methanol  permeation  and  the  pertinent  mixed  potential  effect  in  a  direct  methanol  fuel  cell  (DMFC)  is 
presented.  In  this  model  a  DMFC  is  divided  into  seven  compartments  namely  the  anodic  flow  channel,  the  anodic  diffusion  layer,  the  anodic 
catalyst  layer,  the  proton  exchange  membrane  (PEM),  the  cathodic  catalyst  layer,  the  cathodic  diffusion  layer  and  the  cathodic  flow  channel.  All 
compartments  are  considered  to  have  finite  thickness,  and  within  every  one  of  them  a  set  of  governing  equations  are  given  to  stipulate  methanol 
transport  and  oxygen  transport.  For  the  flow  channels,  fluid  dynamics,  which  could  substantially  lower  the  local  methanol  concentration  within 
catalyst  layers  is  taken  into  account.  With  the  knowledge  of  local  concentrations  of  the  species,  the  electrochemical  reaction  rates  within  both 
catalyst  layers  can  be  quantified  by  a  kinetic  Tafel  expression.  For  the  anodic  catalyst  layer  the  local  external  current  generated  by  methanol 
oxidation  is  computed;  for  the  cathodic  catalyst  layer,  in  addition  to  the  local  external  current  generated  by  oxygen  reduction,  the  local  internal 
current  as  a  result  of  methanol  permeation  is  also  computed.  With  the  information  of  the  local  internal  current,  the  mixed  potential  effect,  which 
is  responsible  for  adversely  lowering  the  cell  voltage  can  be  analyzed. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  advantages  of  direct  methanol  fuel  cells  (DMFC)  over 
hydrogen  fuel  cells  include  easy  storage  of  the  high  energy  den¬ 
sity  liquid  fuel,  direct  fuel  feeding  without  reforming,  and  low 
operating  temperature.  It  is  therefore  considered  by  many  people 
the  most  promising  alternative  power  source  for  mobile  appli¬ 
cations  and  electric  vehicles. 

As  shown  schematically  in  Fig.  1,  the  structure  a  DMFC 
consists  of  seven  major  compartments  namely  the  anodic  flow 
channel,  the  anodic  diffusion  layer,  the  anodic  catalyst  layer, 
the  proton  exchange  membrane  (PEM),  the  cathodic  catalyst 
layer,  the  cathodic  diffusion  layer  and  the  cathodic  flow  chan¬ 
nel.  The  anodic  flow  channel  is  the  passage  of  low  concentration 
methanol  solution.  As  the  solution  is  pumped  through  the  chan¬ 
nel,  a  small  fraction  of  methanol  diffuses  through  the  anodic 
diffusion  layer  and  reaches  the  anodic  catalyst  layer.  Within 


*  Corresponding  author.  Tel.:  +886  3  5742865;  fax:  +886  3  5713849. 
E-mail  address:  chih_hao_chen@yahoo.com  (C.-H.  Chen). 

0378-7753/$  -  see  front  matter  ©  2006  Elsevier  B.V.  All  rights  reserved. 
doi:10.1016/j.jpowsour.2006.03.005 


this  layer,  where  Pt-Ru  is  the  most  widely  used  catalyst  today, 
methanol  oxidizes  and  produces  carbon  dioxide  (CO2)  via  the 
following  reaction: 

CH3OH  +  H20  -*  C02  +  6H+  +6e“. 

CO2  then  diffuses  back  into  the  anodic  flow  channel  and  exits 
with  the  solution.  The  protons,  which  travel  through  the  PEM, 
and  the  electrons,  which  travel  through  some  external  load,  reach 
the  cathodic  catalyst  layer,  where  Pt  is  the  catalyst,  to  undergo 
the  following  half-cell  reaction  with  oxygen  that  comes  from  the 
cathodic  flow  channel: 

3 /202  +6H+  +6e“  —*■  3H20. 

Because  these  electrons  travel  through  the  external  load,  the 
pertinent  current  is  referred  to  as  the  ‘external’  current  in  this 
paper.  The  overall  reaction  can  therefore  be  written  as: 

CH3OH  +  3/202^  C02+2H20. 

Despite  its  advantages  over  hydrogen  fuel  cells,  a  few  engi¬ 
neering  obstacles  of  the  DMFC  remain  to  be  overcome.  On  the 
one  hand,  the  sluggish  catalytic  activity  of  the  anode  makes  a 
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Nomenclature 

Aceii  cell  area 

^ch-Joh  methanol  feed  concentration 

^ch3oh  reference  methanol  concentration 
Cch3oh  local  methanol  concentration 
^ch3oh  gase°us  methanol  concentration  at  saturated 
vapor  pressure 

Ch2o  local  water  concentration 

Cq2  d  oxygen  feed  concentration 

CqI  reference  oxygen  concentration 

Co2  local  oxygen  concentration 

daf  width  of  the  anodic  flow  channel 

dc f  width  of  the  anodic  flow  channel 

A:h3oh,h2o  bulk  diffusion  coefficient  of  methanol 
in  water 

DcH3OH,air  bulk  diffusion  coefficient  of  gaseous 
methanol  in  air 

Z)ch3 oh, pem  diffusion  coefficient  of  methanol  in  PEM 
^o2,air  bulk  diffusion  coefficient  of  oxygen  in  air 

^ch3oh  Ho  o  effective  diffusion  coefficient  of  methanol  in 
diffusion  layers 

^ch3  oh  Ho  o  effective  diffusion  coefficient  of  methanol  in 
the  anodic  catalyst  layer 

^CH3OH  Ho  o  effective  diffusion  coefficient  of  methanol  in 
the  cathodic  catalyst  layer 

d  ^jjf 

Dq>1  air  effective  diffusion  coefficient  of  oxygen  in  diffu¬ 
sion  layers 

D0i  air  effective  diffusion  coefficient  of  oxygen  in  the 
cathodic  catalyst  layer 

E  difference  of  electrode  potentials  in  a  DMFC 
/af  flow  rate  of  anodic  flow  channel 
/cf  flow  rate  of  cathodic  flow  channel 
ycee„  cell  external  current  density 

j^eii  cell  internal  current  density 

F  Faraday’s  constant 

jo2  local  current  density  from  oxygen  reduction 
Jc h3oh  l°cal  external  current  density  from  methanol  oxi¬ 
dization 

cL  ref 

j0’Ch3oh  reference  exchange  current  density  of  methanol 
in  the  anode 

j0’ch3oh  reference  exchange  current  density  of  methanol 
in  the  cathode 

j'q  local  external  current  density  from  oxygen  reduc¬ 

tion 

j0’q2  reference  exchange  current  density  of  oxygen  in 
the  cathode 

cL  ref 

/c’h3oh  reference  catalyst  layer  thickness  for  methanol  in 
the  anode 

c  ref 

/c’h3oh  reference  catalyst  layer  thickness  for  methanol  in 
the  cathode 

Iq  reference  catalyst  layer  thickness  for  oxygen  in 
the  cathode 

/af  thickness  of  the  anodic  flow  channel 


/ad  thickness  of  the  anodic  diffusion  layer 
/ac  thickness  of  the  anodic  catalyst  layer 
/m  thickness  of  the  PEM 
/cc  thickness  of  the  cathodic  catalyst  layer 
led  thickness  of  the  cathodic  diffusion  layer 
lcf  thickness  of  the  cathodic  flow  channel 
Mch3oh  molecular  weight  of  methanol 
Mh2o  molecular  weight  of  water 
^ch3oh  number  of  transferred  electrons  per  methanol 
molecule 

no2  number  of  transferred  electrons  per  water 
molecule 

Nch3oh  local  methanol  flux 
No2  local  oxygen  flux 

Nu2o  local  water  flux 

rceii  cell  interfacial  resistence 

Vceii  cell  output  voltage 

^CH3OH,a  decay  width  of  methanol  concentration  along 
the  anodic  flow  channel 

^ch3oh,c  decay  width  of  methanol  concentration  along 
the  cathodic  flow  channel 

wq2  c  decay  width  of  oxygen  concentration  along  the 
cathodic  flow  channel 

Greek  letters 

a\  CH3oh  an°dic  transfer  coefficient  of  methanol  in  the 
anode 

at  ch3oh  an°dic  transfer  coefficient  of  methanol  in  the 
cathode 

ac  ch3  oh  cathodic  transfer  coefficient  of  methanol  in  the 
anode 

ac  ch3  oh  cathodic  transfer  coefficient  of  methanol  in  the 
cathode 

aca  q2  anodic  transfer  coefficient  of  oxygen  in  the  cath¬ 
ode 

acc  0o  cathodic  transfer  coefficient  of  oxygen  in  the  cath¬ 

ode 

£d  void  fraction  of  diffusion  layers 
£^c  volume  fraction  of  solid  phase  in  anodic  catalyst 
layer 

£gC  volume  fraction  of  solid  phase  in  cathodic  catalyst 

layer 

volume  fraction  of  ionomer  phase  in  catalyst  lay¬ 
ers 

yo2  reaction  order  of  oxygen 
^ch3oh  reacti°n  order  of  methanol  in  the  anode 
kcH3OH  reacti°n  order  of  methanol  in  the  cathode 
Yq2  reaction  order  of  oxygen  in  the  cathode 
^ch3oh  l°cal  external  overpotential  of  methanol 
rj o2  local  overpotential  of  oxygen 
rf0  local  external  overpotential  of  oxygen 
Km  specific  conductivity  of  the  PEM 
^.ac,eff  effective  conductivity  of  solid  phase  in  the  anodic 
catalyst  layer 
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,xc,eff 

acs 

effective  conductivity  of  solid  phase  in  the 
cathodic  catalyst  layer 

^c,eff 

effective  conductivity  of  ionomer  phase  in  cata¬ 
lyst  layers 

<c 

conductivity  of  solid  phase  in  the  anodic  catalyst 
layer 

<c 

conductivity  of  solid  phase  in  the  cathodic  catalyst 
layer 

conductivity  of  ionomer  phase  in  catalyst  layers 

^h2o 

drag  coefficient  of  water 

PCH30H 

density  of  methanol 

Ph2o 

density  of  water 

higher  methanol  concentration  favorable.  On  the  other  hand  the 
methanol  permeation  problem,  which  does  not  exist  in  hydro¬ 
gen  fuel  cells,  generates  a  mixed  potential  at  the  cathode  and 
adversely  lowers  the  cell  output  voltage  at  higher  methanol 
concentration.  Carbon  monoxide  (CO)  poisoning  poses  another 
problem  by  reducing  the  electrocatalytic  activity  and  active  area 
of  the  catalyst.  Another  important  issue  is  that  methanol  trans¬ 
port  may  be  hindered  by  CO2  that  diffuses  back  into  the  anodic 
flow  channel  after  being  released  within  the  anodic  catalyst 
layer. 

Methanol  permeation  is  considered  the  most  serious  prob¬ 
lem.  The  PEM  is  fully  or  almost  fully  hydrated  when  a  DMFC 
is  operating  and  methanol  can  therefore  diffuse  through  the  PEM 
and  reach  the  cathodic  catalyst  layer.  The  catalyst  in  the  cathode, 
which  is  made  of  Pt  and  adopted  to  reduce  oxygen,  also  oxidizes 
the  permeating  methanol  at  the  same  time.  Since  the  exchanged 
electrons  in  these  redox  reactions  do  not  go  through  any  exter¬ 
nal  load,  we  refer  to  this  current  as  the  ‘internal’  current.  How 
this  internal  current  influences  the  performance  of  a  DMFC  is 
explained  in  the  following.  The  cell  output  voltage  V  is  related 


to  the  potential  difference  E  by: 

V  =  E  -  r)a-  rjc  ~  IR, 

where  r /a  and  rjc  are  the  activation  overpotentials  of  the  anode 
and  the  cathode,  respectively,  and  IR  is  the  ohmic  loss.  The 
cathodic  catalyst  layer  is  assumed  to  be  of  zero  thickness.  Since 
the  total  amount  of  reduced  oxygen  must  account  for  both  the 
internal  current  density  i1  and  the  external  current  density  /e,  r]c 
is  stipulated  by  the  following  Buttler-Volmer  equation: 

Oc 


where  ba  and  bc  are  the  Tafel  slopes  of  oxygen  oxidation  and 
reduction,  respectively.  If  we  keep  f  constant  and  raise  i\  rjc  will 
be  raised  and  therefore  V  will  be  lowered.  Because  contributions 
from  i1  and  f  to  r]c  are  mixed  together  and  cannot  be  separated, 
77c  is  usually  referred  to  as  the  mixed  potential. 

Performance  of  a  DMFC  relies  on  a  vast  number  of  param¬ 
eters,  including  the  methanol  feed  concentration,  efficiencies 
of  methanol  transport  and  oxygen  transport  within  the  compart¬ 
ments,  the  release  rate  of  gaseous  CO2  and  its  effect  on  methanol 
transport,  the  specific  area  of  catalyst  in  the  catalyst  layers,  the 
thickness  of  the  compartments,  the  impedance  of  the  catalyst 
layer,  the  impedance  of  the  membrane,  the  design  of  flow  chan¬ 
nels,  the  rate  of  methanol  permeation  and  so  on.  Investigating 
the  impact  of  these  parameters  one  by  one  through  experiments 
is  not  time  or  cost  efficient.  In  order  to  help  understand  the 
operation  of  a  DMFC  and  locate  the  key  parameters  on  cell  per¬ 
formance,  a  theoretical  model  is  essential.  A  number  of  such 
models  already  exist  in  the  literature. 

Scott  et  al.  [1]  presented  a  DMFC  model  which  accounts  for 
the  influence  of  methanol  transport  and  of  CO2  flow.  The  model 
considers  hydrodynamics  in  the  flow  channels  associated  with 
the  production  of  CO2  (in  an  empirical  way)  and  methanol  trans¬ 
port  within  the  anodic  diffusion  layer,  the  anodic  catalyst  layer 
and  the  PEM.  Within  the  PEM  water  transport  is  accounted  for 
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Fig.  1.  Schematic  of  the  DMFC  which  is  divided  into  seven  compartments  namely  the  anodic  flow  channel,  the  anodic  diffusion  layer,  the  anodic  catalyst  layer,  the 
PEM,  the  cathodic  catalyst  layer,  the  cathodic  diffusion  layer  and  the  cathodic  flow  channel. 
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by  diffusion  and  electro-osmotic  drag.  Within  the  anodic  catalyst 
layer,  methanol  oxidation  is  accounted  for  by  Tafel  type  kinet¬ 
ics.  Experimental  data  were  used  for  validation.  Mass  transport 
in  the  anode  is  found  to  be  a  factor  that  may  limit  the  perfor¬ 
mance  of  the  fuel  cell.  The  open  circuit  potential  as  a  result 
of  methanol  permeation  is  incorporated  with  a  semi-empirical 
equation  instead  of  theoretical  analysis. 

Jeng  and  Chen  [2]  reported  a  model  for  the  anode  of  a  DMFC. 
In  the  anodic  flow  channel,  methanol  transport  associated  with 
the  production  of  CO2  is  accounted  for  in  an  empirical  way. 
Diffusion  and  convection  account  for  methanol  transport  and 
water  transport  in  the  diffusion  layer  as  well  as  in  the  cata¬ 
lyst  layer.  Within  the  PEM  water  transport  is  accounted  for  by 
diffusion  and  electro-osmotic  drag  while  methanol  transport  is 
accounted  for  in  a  similar  way  as  in  the  diffusion  layer.  The 
electrochemical  reaction  in  the  catalyst  layer  is  stipulated  by 
a  kinetic  Tafel  expression.  The  mixed  potential  effect  is  not 
addressed. 

Meyers  and  Newman  [3]  developed  a  thermodynamic  frame¬ 
work  with  which  they  are  able  to  model  the  multi-component 
transport  of  species  within  the  PEM  [4] .  In  this  model  the  authors 
quantify  the  gradients  of  electrochemical  potential  to  describe 
the  driving  forces  for  the  multi-component  transport.  In  addition, 
a  kinetic  model  considering  multiple  reaction  steps  is  developed 
to  describe  methanol  oxidation  on  the  catalysts.  With  this  model, 
a  simulation  of  the  direct  methanol  fuel  cell  is  carried  out  and 
general  aspects  of  its  design  are  quantified  [5].  The  mixed  poten¬ 
tial  effect  is  again  not  addressed. 

In  Divisek  et  al.  [6],  the  authors  developed  and  tested  a  two- 
dimensional  model  of  a  DMFC.  In  their  model,  the  water  trans¬ 
port  is  assumed  to  follow  a  standard  two-phase  flow  mechanism, 
which  can  be  seen  as  a  generalized  Darcy’s  law.  The  interac¬ 
tions  between  the  gases  and  the  ones  between  the  gases  and  the 
pore  walls  are  also  taken  into  consideration.  The  mass  transport 
of  dissolved  species  is  accounted  for  with  a  standard  transport 
model.  For  charged  species,  appropriate  potential  equations  are 
used.  Energy  considerations  are  derived  from  the  Fourier  law, 
with  convection  caused  by  the  flow  of  fluids  and  gases.  Electro¬ 
chemical  reactions  are  accounted  for  with  methanol  and  oxygen 
kinetics.  Condensation  and  evaporation  processes  are  also  con¬ 
sidered.  Although  methanol  permeation  is  considered,  the  mixed 
potential  effect  is  not  addressed. 

In  Wang  and  Wang  [17],  two-phase  mass  transfer  is  con¬ 
sidered  in  the  cathode  (water  and  oxygen)  and  in  the  anode 
(methanol  solution  and  CO2)  of  a  DMFC.  Oxygen  reduction  and 
methanol  oxidation  are  accounted  for  by  Tafel  kinetic  equations. 
Simulation  results  show  that  cell  performance  relies  heavily  on 
methanol  feed  concentrations.  With  small  methanol  feed  con¬ 
centrations  the  cell  suffers  from  low  limiting  current  densities. 
For  feed  concentrations  below  1  M,  an  increase  in  methanol  con¬ 
centration  leads  to  a  small  decrease  in  cell  voltage.  For  a  feed 
concentration  larger  than  2  M,  the  cell  voltage  is  reduced  sub¬ 
stantially  by  excessive  methanol  permeation  and  the  maximum 
current  is  limited  by  the  oxygen  depletion  at  the  cathode.  The 
mixed  potential  effect  as  a  result  of  methanol  permeation  is  taken 
into  account  by  assuming  that  the  methanol  transported  to  the 
cathode  is  completely  depleted. 


The  mixed  potential  effect  is  not  easy  to  simulate  and,  in 
the  models  mentioned  above,  is  unaddressed,  calculated  in  an 
empirical  way  or  handled  with  a  simple  assumption  of  full 
depletion  of  the  permeating  methanol  by  the  cathode.  Empirical 
approaches  are  often  useful  in  correlating  experimental  data  if 
the  model  contains  sufficient  insights  of  the  system,  but  are  less 
helpful  on  the  investigation  of  cell  parameters  or  on  the  effects 
of  changing  cell  designs.  The  assumption  of  full  depletion  of 
the  permeating  methanol  by  the  cathode  may  not  always  be 
true,  especially  when  the  catalytic  activity  of  the  cathode  is 
not  so  efficient,  as  found  out  and  discussed  in  the  Results  and 
Discussion  section.  The  mathematical  model,  which  is  proposed 
here  is  based  upon  the  description  of  the  physicochemical  pro¬ 
cesses  which  dictate  the  behavior  of  electrochemical  systems, 
namely,  mass  transport  and  reaction  kinetics.  The  features  that 
distinguish  our  model  from  the  foregoing  ones  are: 

1 .  fluid  dynamics  of  the  flow  channels  which  may  influence  the 
concentrations  of  species  at  the  catalyst  layers,  change  the 
reaction  rates  and  impact  the  performance  of  the  cell, 

2.  mass  transport  and  reaction  kinetics  of  the  cathode,  and 

3.  a  means  to  estimating  the  intensity  of  the  internal  current,  the 
mixed  potential  effect  and  consequently  the  impact  on  cell 
output  voltage. 

2.  Theory  and  calculation 

Our  model  is  based  on  the  following  set  of  assumptions: 

1.  The  fuel  cell  is  operated  isothermally  at  30° C  in  a  steady 
state. 

2.  There  is  no  pressure  difference  between  the  compartments. 

3.  Methanol  flux  into  the  anodic  flow  channel  is  much  greater 
than  methanol  flux  into  the  anodic  diffusion  layer.  This 
ensures  that  methanol  concentration  variation  is  small  along 
the  channel.  The  same  assumption  applies  to  oxygen  flux  in 
the  cathodic  flow  channel. 

4.  The  effects  of  generated  products,  carbon  dioxide  and  water, 
on  methanol  transport  and  on  oxygen  transport  are  neglected. 

Furthermore,  as  will  be  discussed  in  the  next  section,  the 
reduction  in  the  electrocatalytic  activity  and  active  reaction 
area  of  the  catalyst  due  to  CO  poisoning  is  accounted  for 
by  selecting  appropriate  reference  exchange  current  densities, 

namely 7q’^oh,  Jo,ch3OH  and  kjf  However,  it  should  be 
noted  that,  since  we  only  discuss  steady- state  operations,  the 
transient  effect  of  CO  poisoning  on  these  exchange  current  den¬ 
sities  is  not  included  in  this  model. 

A  variable  map  is  presented  in  Table  1 .  The  concentration  and 
flux  of  methanol,  Cch3oh  and  AAh3oh,  are  considered  all  over 
the  cell.  The  concentration  and  flux  of  oxygen,  Cq2  and  No2 ,  are 
considered  only  in  the  cathode  because  oxygen  does  not  pene¬ 
trate  the  PEM.  Within  the  anodic  catalyst  layer  where  methanol 
oxidizes,  the  external  current  density  from  methanol  oxidation 
7ch3oh  and  its  activation  overpotential  0ch3ou  are  added. 
Within  the  cathodic  catalyst  layer,  in  addition  to  the  external 
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Table  1 

Variable  map  showing  the  variables  which  are  solved  for  in  the  seven  compartments  of  the  DMFC 
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V 
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V 

4 

V 

jo2 

V 

current  density  from  oxygen  reduction  Jq  and  its  activation 
overpotential  rfo2 » the  total  current  density  from  oxygen  reduc¬ 
tion  jo2  is  needed  to  account  for  the  internal  currents.  Our  choice 
differs  from  other  approaches  in  the  literature  and  substantially 
simplifies  the  governing  equations  and  the  boundary  conditions. 

2.7.  Governing  equations 

The  following  notations  are  used  to  simplify  the  govern¬ 
ing  equations.  7?S(C ,  q)  is  the  abbreviation  of  the  kinetic  Tafel 
expression: 


where 


AfConv 

7VCH3OH 


and 


Mh2oCch3OH  \x=xad  Nr2q  I x=x ad 
Ph2o 


„diff  r,  dCcil.OII 

^CHjOH  =  -ArH30H,H20 - — - • 

The  last  equation  can  be  reorganized  as: 
dCcH3OH 


Afdiff 

7VCH3OH 


dx 


A:h3oh,h2o 


BUC,  rj)  = 


■l,ref 

Jq,$ 

,l,ref 

*S 


c 


rl 


^bref 

Vs 


X 


exp 


al,nsF 


RT 


q 


exp 


,1 

'C,S* 


ac  F 


RT 


where  s  stands  for  species  (CH3OH  or  O2)  and  1  for  location  (a 
for  anode  or  c  for  cathode).  Effective  diffusion  coefficients  (Ds) 
and  effective  conductivities  (ks)  are  defined  as: 


This  reorganization  is  also  applied  to  mass  transfer  in  all  the 
other  compartments. 

By  fluid  dynamics  and  mass  transfer  calculation  in  Appendix 
A,  0H  must  satisfy  the  following  equation: 


d/vdiff 

q7VCH3OH 

dx 


6/af 


o 

^CH3OH,a^af4f 


(X  -  Xaf)(/af  -  X  +  Xaf)CcHsOH 


(2) 


r-vd,eff  _  /  dA/2 

^CH30H,H20  —  w  )  7>>ch3OH,H20, 

nac’eff  _  oac  \ 3/2  n 

^CH30H,H20  —  G  Ss  £m )  77ch30H,H20, 

7^cc,cff  /  \  cc  c  \3/2 

^CH30H,H20  =  G  —  £s  —  £m>  f7cH30H,H20 

AC2eA,r  =  (|-£sC-4)3/2Oo2.A,r. 


ac,eff  _  /  ac\3/2  ac 
s  ■ 


4f  =  (4)3/2 


m 


K 


m  5 


<c’eff  =  (ef  )3/Vc 


S  ’ 


using  Bruggeman’s  correction  [7-9].  The  derivation  of  the  gov¬ 
erning  equations  in  the  seven  compartments  of  a  DMFC  is  shown 
below. 


2.7.2.  The  anodic  diffusion  layer 

In  this  compartment,  methanol  transfer  is  driven  by  diffusion 
and  convection: 

dCcH3OH  _  Nch3oh  ^7h2o^h2o^ch3oh 

TZ  _  nd,eff  1  TZeff  '  ^ 

^ch3oh,h2o  ^ch3oh,h2o^h2o 

Because  no  electrochemical  reactions  occur,  the  methanol  flux 
and  the  water  flux  are  constant: 


dA^CH3OH 

dx 


dNu2o  _  Q 
dx 


2.7.7.  The  anodic  flow  channel 

In  this  compartment  methanol  flux  results  from  diffusion  and 
convection: 

AT  —  Afdiff  1  A/-COnV 

IXCU3OU  —  NcH3OH  +  ^CH3OH’ 


2.1.3.  The  anodic  catalyst  layer 

In  this  compartment,  methanol  transfer  is  a  result  of  diffusion 
and  convection: 


dCcH3OH 

dx 


Ach3oh  ^H2oAfit2oCcH3OH 

T^ac^eff  ^ac,eff 

^ch3  oh,  h2o  ^ch3  oh,  h2o  ^h2o 
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Methanol  oxidation  which,  can  be  stipulated  by  a  kinetic  Tafel 
expression  changes  the  methanol  flux  and  the  water  flux: 


djcH3OH 

dx 

-  j6ch3oh((7ch3oh,  ^?CH3Oh)’ 

(7) 

dNcH3OH 

_  1  djcH3OH 

(8) 

dx 

ncn3ouE  dx 

dA^H2o  

dNcu3ou 

(9) 

dx 

dx 

^CH3OH’  external  overpotential  by  methanol,  is  changed  by 
the  electron  flow  in  solid  phase  and  the  proton  flow  in  ionomer 
phase  [2]: 

^^ch3oh  f  1 
dx  ~~ 


+  j  7ch3oh 


1 


ac,eff 


7-e 

J  cell ' 


(10) 


2.1.4.  ThePEM 

In  this  compartment,  methanol  transfer  is  caused  by  diffusion 
and  convection: 


dCcH3OH 

dx 


Nch3oh 

£>ch3oh,pem 


+ 


Mh2  O  Nh2  O  c  ch3  oh 
£>ch3oh,pemPh2o 


(11) 


No  electrochemical  reaction  occurs  in  this  compartment  to 
change  the  methanol  flux  or  the  water  flux: 


dNcH3OH 

dx 


dA^H2o 

dx 


(12) 

(13) 


We  further  assume  that  the  PEM  is  fully  hydrated  so  that  water 
concentration  is  constant  and  that  no  water  diffusion  takes  place. 


dA^CH3OH 

1  /  d/o2 

d4  \ 

(19) 

dx 

^ch3ohE  V  dv 

dx  J  ’ 

dA4t2o  

dNcH3OH  0dA^o2 

J 

(20) 

dx 

dx  dx 

Once  again  the  electron  flow  in  the  solid  phase  and  the  proton 
flow  in  the  ionomer  phase  change  the  external  overpotential  by 
oxygen  if(h : 


dl1o2  =  (  1  |  1 

dx  V/^eff  4c’eff 

The  sum  of  the  methanol  internal  current  density  and  the  oxygen 
internal  current  density  must  be  zero.  Therefore,  the  following 
equation  must  hold  all  over  this  compartment: 


d,/o2 

dx 


d 4 

dx 


+  5CH3Oh(CcH3OH,  E  +  T]Q2)  =  0. 


(22) 


2.1.6.  The  cathodic  diffusion  layer 

In  this  compartment,  methanol  transfer  is  due  to  diffusion 
and  convection: 


dCcH3OH 

dx 


D 


Nch3oh  ,  ^H2oA^H2O^CH3OH 

nd,eff 

^ch3oh,h2o^h2o 


d,eff 

ch3oh,h2o 


+ 


Oxygen  transfer  is  a  result  of  diffusion: 


dCp2 

dx 


Nq2 

nd,eff 

^02,air 


(23) 


(24) 


No  electrochemical  reaction  occurs  to  change  the  methanol  flux, 
the  oxygen  flux  or  the  water  flux: 


2.1.5.  The  cathodic  catalyst  layer 

In  this  compartment,  methanol  transfer  is  due  to  diffusion 
and  convection: 


dCcH3OH 

dx 


D 


Nch3oh  ,  ^h2oNh2oCch3oh 

r^cc,eff 

^ch3oh,h2o^h2o 


cc,eff 


+ 


ch3oh,h2o 

Oxygen  transfer  is  a  result  of  diffusion 


dCp2 

dx 


N02 

^cc,eff 

^02,air 


(14) 


(15) 


Oxygen  reduction  generates  both  the  external  current  and  the 
internal  current.  The  combination  of  both  and  the  former  along 
are  stipulated  by  kinetic  Tafel  expressions: 


d,/o2 

dx 


dx 


Bc02(Co2,  tjo2), 

Bq2(Co2,  0q2)- 


(16) 

(17) 


The  methanol  flux,  the  oxygen  flux  and  the  water  flux  all  change 
accordingly  with  the  generation  of  these  current  densities: 


dAfp2  1  djp2 

dx  no2F  dx 


(18) 


dNcH3OH 

dx 


dA^H2o  _  Q 
dx 


(25) 

(26) 
(27) 


2.1.7.  The  cathodic  flow  channel 

In  this  compartment,  methanol  transfer  and  oxygen  transfer 
can  be  derived  as  in  the  anodic  flow  channel: 


dNcH3OH 

dx 


6/cf 
01 CH3  OH ,  c 


(x  xcf)(xend 


x)Och3oh,(28) 


dCcH3OH 

dx 


Nch3oh 

7^CH3OH,air  ’ 


(29) 


dNp2 

dx 


6/cf 

OJQ2,cdcfl^f 


(v  xcf)(xend 


x)Cq2  , 


(30) 


dCo2  _  No2 

dx  D  02,air 


(31) 


Note  that  methanol  is  in  gas  phase  here. 
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2.2.  Boundary  conditions 

where  can  be  determined  by 

A  set  of  boundary  conditions  are  required  to  solve  the  differ¬ 
ential  equations  above. 

47ch3oh<7ch3oh  4^h2o<7h2o  _  i 

!  T  • 

Pch3oh  Ph2o 

2.2.1.  Atx-xaf 

^feed 

EcH3OH  —  ECH3oh> 

(32) 

2.2.8.  Atx  =  xenc[ 

^feed 

co2  —  c0l  , 

No2  =  0, 

(48) 

(49) 

yydiff  _  a 

7Vch3oh  —  u- 

(33) 

6"ch3oh  =  0, 

(50) 

2.2.2.  Atx-xa(i 

Nch3oh  =  0. 

(51) 

6"ch3oh  and  AAh3oh  arecontinuous. 

(34) 

2.3.  Numerical  method 

Water  flux  results  from  water  consumption  by  the  anodic  catalyst 
layer  and  electro-osmosis  in  the  PEM: 

Eqs.  (1)— (3 1)  form  a  set  of  first-order  differential 

equa- 

Nh2o  = 


1 


/lCH3OH 


+  ^h2o 


7e 

J  cell 

F 


(35) 


2.2.3.  Atx-xac 

Cch3oh  and  A^ch3oh  arecontinuous, 


Jch.oh 


=  0. 


(36) 

(37) 


2.2.4.  At  x-  xm 

Cch3oh  and  A^ch3oh  arecontinuous, 


7CH3OH  —  ./cell' 


2.2.5.  Atx-xcc 

6"ch3oh  and  AAh3oh  arecontinuous, 
^Vo2  =  0, 
jo2  =  0. 
fo2  =  0. 


(38) 

(39) 


(40) 

(41) 

(42) 

(43) 


Eq.  (41)  comes  from  our  assumption  that  oxygen  does  not  pen¬ 
etrate  the  PEM  and  exists  only  in  the  cathode. 

2.2.6.  Atx-xC(i 

Cch3oh,  Nch3oh,  Cq2  and  Nq2  arecontinuous,  (44) 


•e  _  -e 

J  02  —  “7cell 


(45) 


2.2.7.  Atx-xcf 

Nch3oh,  Cq2  and  Nq>2  arecontinuous. 


(46) 


Methanol  vaporizes  at  this  boundary  and  its  concentration  in  gas 
phase  can  be  determined  by  its  molar  fraction  in  the  solution: 


<7ch3oh 

6"ch3oh 

^vap 

_cch3oh_ 

X>Xcf 

_<7ch3oh  +  Ch2o. 

(47) 


X<Xcf 


tions,  which  can  be  readily  solved  with  boundary  conditions  Eq. 
(32)— (5 1)  by  applying  Runge-Kutta  method  of  order  4.  Because 
the  boundary  conditions  are  not  at  the  same  location,  the  solv¬ 
ing  procedure  must  include  first  guesses  of  variables  and  a  large 
number  of  iterations  to  make  the  solution  converge. 

2.4.  Determination  of  cell  voltage 
The  total  power  density  of  the  cell  is 

Pi  ^O'cell  3”  7cell)* 

The  usable  power  density  of  the  cell  is 
Pn  =  V/cell  ■ 

The  power  density  of  anodic  activation  is 
_  fXm  e  d/CH,OH 

P  aa  —  /  *7cH3OH  ^  ax- 

^  X3l C 

The  power  density  of  cathodic  activation  is 


Pac=  /  (E  +  rjo2) 


d/o2  djo2  V..  ,  fXm  ..  dio2 


dx 


dx 


dx  + 


oo2 


dx 


dx, 


Xac  \  /  J  Xac 

where  the  first  term  is  for  methanol  and  the  second  term  for 
oxygen.  The  power  density  of  IR  loss  is 


PlR  = 


Och3oh 


x 


ac 


c,eff 


^  _j_  ^cell  7ch3Oh)^ 


ac,eff 


,JCcd  l  O'cell  +  fQf  ,  (i02)2 


dx 


+ 


lcc 


(7’cell)2 


c,eff 

Km 


+ 


o2 


cc,eff 

Ks 


dx 


+ 


K 


lm  (jCell)  ^celRcell’ 


m 


where  the  first  term  is  for  the  anodic  catalyst  layer,  the  second 
term  for  the  cathodic  catalyst  layer[2] ,  the  third  term  for  the  PEM 
and  the  forth  term  for  interfacial  resistance.  By  conservation  of 
power  density, 

Pu  =  Pt~  Paa  -  Pac  “  PlR, 
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we  then  have 


^cellj'cell  £(Jccll  ./cell)’ 


m 


^CthOH 


^^CH3OH 


-ac 


dx 


fX  m  / 

/  (F  +  7?02)  ( 

Xrc  \ 


d./'o2  dy0. 


X 


lac 

*-^m 

ac 


dx 


dx 


dx 


m2 


d./o2 

dx 


dx, 


‘X 


m 


0ch3OH 


•^ac 


c,ef  / 
/Cm 


^  _l_  ^ceii  jch3oh)2 


dx, 


ac,eff 
K  s 


dx 


^cd  /  ( 7®ell  +  j‘o9 )  Ooo ) 


•*cc 


c,eff 


)2'  _|_  kJo2 


cc,eff 

Ks 


dx, 


Ocell)” 


K 


Ocell)  ^cell^cell* 


m 


After  rearranging  the  terms  and  applying 


’■^m  17  I  "J02 


L  E 

u  Aac 


<U? 


dx 


djO; 

dx 


dx  =  E j*ell  we  have 


^cellj'cell  ^v/cell’ 


’•Wn 


•^ac 


e  dy'cH,0H  , 
'tiiion  ^  dx- 


‘X 


m 


djf, 

m2  ~ r-1  dx> 


Xac 

>X™  /  ^  ’e 


dx 


Lm 


0 


CH3Oh)~  Ocell  iCH3OH)2 


x 


ac 


c,eff 


ac,eff 

Ks 


dx 


^cd  (  (j®ell  +  j‘q2)2  ^  (jp2)2  |  ^ 


XCc 


c,eff 


cc,eff 


Ocell)2 


/C 


/m  Ocell )  ^cell^cell' 


m 


(52) 


3.  Results  and  discussion 

In  order  to  validate  our  model,  a  cell  is  assembled  in  our  labo¬ 
ratory  with  which  experimental  data  is  generated  and  compared 
to  model  predictions.  The  catalyst  layers  of  the  cell  are  made  of 
Nation  117,  Teflon-treated  carbon  clothes  and  commercial  cata¬ 
lysts  from  Johnson  Matthey  with  these  weight  ratios:  Pt  20  wt.% 
and  Ru  10  wt.%  for  the  anode,  Pt  20%  for  the  cathode.  Amounts 
of  Pt  usage  on  catalyst  layers  are  0.9  mg  cm-2  for  both  the 
anode  and  the  cathode.  A  membrane  made  of  Nation  1 17  is  hot 
pressed  between  the  catalyst  layers  at  130  °C  to  form  the  mem¬ 
brane  electrode  assembly  (ME  A).  The  ME  A  is  then  sandwiched 
between  two  graphite  plates  on  which  flow  channels  are  carved. 
After  the  cell  is  assembled,  we  delineate  its  I-V  characteris¬ 
tic  curves  at  different  methanol  feed  concentrations  as  shown 
in  Fig.  2. 

Next  we  determine  key  parameter  values  of  the  model  through 
calibration,  by  fitting  model-predicted  I-V  characteristic  curves 


Fig.  2.  Polarization  curves  of  our  cell  compared  with  model  predictions  after 
calibration. 

to  experimental  data.  The  best  fitting,  also  shown  in  Fig.  2,  is 
obtained  with  the  set  of  parameter  values  listed  in  Table  2. 

With  key  parameter  values  determined,  the  mixed  potential 
effect  can  be  analyzed  by  examining  the  difference  between  the 
cathodic  activation  overpotentials  of  oxygen  with  the  internal 
current  enabled  and  disabled.  In  the  case  of  disabled  internal 
current,  the  way  methanol  transport  to  the  cathode  is  not  changed 
while  methanol  oxidation  within  the  cathodic  catalyst  layer  is 
assumed  not  to  take  place.  This  difference,  as  shown  in  Fig.  3 
where  the  methanol  feed  concentration  is  still  1  M,  decreases 
as  cell  current  density  increases.  This  can  be  understood  since 
higher  cell  current  densities  gives  lower  methanol  flux  into  the 
cathode,  as  shown  in  Fig.  5,  which  will  be  discussed  later.  This 
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Fig.  3.  Comparison  of  mixed  potentials  with  internal  current  enabled  and  dis- 
abled  at  Cg£0H  =  1  M. 
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Table  2 

Parameter  values  used  in  numerical  simulations 


Parameter 

Value 

Reference 

Aceii  (cm2) 

4 

Measurement 

A2h3oh,H20  (cm2  s-1) 

1.93  x  10“5 

[10] 

£>CH3OH,air  (CUlV1) 

1.569  x  10-> 

[11] 

A:h3oh,pem  (cm2  s-1) 

4.9  x  10"6 

[13] 

£>02,air  (cm2  S_1) 

1.02  x  101 

[10] 

daf  (cm) 

0.1 

Measurement 

dcf  (cm) 

0.1 

Measurement 

Jaf  (cm3  s_ 1 ) 

0.83 

Measurement 

Jcf  (cm3  s—  1 ) 

1.66 

Measurement 

4;ch3oh  (A cm-2) 

4.5  x  10“4 

Calibration 

ioiclon  (A  cm"2) 

4.5  x  10-4 

Calibration 

'ch3oh  (cm) 

0.03 

Measurement 

'ch'oh  (cm) 

0.03 

Measurement 

foZ  (A  cm-2) 

1  x  10“4 

Calibration 

lC0f  (cm) 

0.03 

Measurement 

/af  (cm) 

0.1 

Measurement 

/ad  (cm) 

0.03 

Measurement 

kc  (cm) 

0.001 

Measurement 

/m  (cm) 

0.015 

Measurement 

kc  (cm) 

0.001 

Measurement 

kd  (cm) 

0.03 

Measurement 

kf  (cm) 

0.1 

Measurement 

4n  (°hm) 

1.4 

Calibration 

aa,CH3OH 

0.153 

Calibration 

Q/C 

ua,CH3OH 

0.153 

Calibration 

uc,CH3OH 

0.12 

Calibration 

cP 

uc,CH3OH 

0.12 

Calibration 

aa,02 

0.0669 

Assumption 

ac,02 

0.0669 

[14] 

^CH3OH 

1.4 

Calibration 

^ch3oh 

2 

Calibration 

ro2 

1 

Calibration 

k^c  (S  cm  1 ) 

8.13  x  106 

[15] 

k “  (S  cm- 1 ) 

8.13  x  106 

Assumption 

Kcm  (S  cm-1) 

1.416  x  KT1 

[13] 

Km  (S  cm-1) 

8.3  x  1(T2 

[16] 

£d 

7.06  x  10-1 

[2] 

Pac 

6  x  10"1 

[2] 

pCC 

6  x  10_1 

Assumption 

£C 

cm 

8  x  10“2 

[2] 

^H20 

2.36 

[12] 

also  indicates  that  methanol  permeation  has  a  weaker  impact  on 
cell  output  voltage  at  higher  cell  current  densities. 

The  model  can  further  be  used  to  predict  methanol  concen¬ 
tration  distributions  and  methanol  flux  inside  the  cell  at  different 
cell  current  densities,  as  shown  in  Figs.  4  and  5  where  the 
methanol  feed  concentration  is  fixed  at  1  M.  Because  it  is  dif¬ 
ficult  to  see  these  curves  within  the  very  thin  cathodic  catalyst 
layer,  we  show  the  magnified  versions  in  the  lower  frames.  Fig.  4 
shows  that  the  methanol  concentrations  fall  substantially  in  the 
anodic  flow  channel  just  before  the  anodic  diffusion  layer.  This 
result,  which  is  just  as  expected,  justifies  the  integration  of  flow 
channel  calculation  into  our  model.  Regarding  the  cathodic  flow 
channel,  due  to  the  high  diffusion  coefficient  of  oxygen  in  air, 
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Fig.  4.  Model-predicted  methanol  concentration  distributions  at  different  cell 
current  densities  with  C§^0H  =  1  M. 


the  model  predicts  little  change  in  oxygen  concentration.  There¬ 
fore,  we  do  not  present  oxygen  concentration  distributions  here. 
Because  electrochemical  reactions  only  occur  in  the  catalyst  lay¬ 
ers,  methanol  flux  is  constant  in  the  diffusion  layers  and  the  PEM, 
as  shown  in  Fig.  5.  We  also  note  that  higher  cell  current  densi¬ 
ties  give  lower  methanol  permeation  rates  and  thus  better  fuel 
efficiencies. 

For  this  particular  cell,  the  catalyst  efficiency  on  methanol 
oxidation  is  considered  poor  but  still  good  enough  to  deplete 
most  of  the  methanol  that  permeates  to  the  cathode,  as  shown  in 
Fig.  5.  However,  the  fraction  of  methanol  that  escapes  into  the  air 
would  increases  if  the  catalyst  efficiency  is  reduced.  To  manifest 
this  phenomenon  with  our  model,  we  deliberately  change  the 
catalyst  efficiency  of  the  cathode  by  defining  the  following  cases: 

CASE_1  Jo  ch3oh  anc^  Jo  02  are  as  listed  in  Table  2. 
CASE_0.25  Jo  CH3OH  anc^  Jo  02  are  multiplied  by  0.25. 
CASE_4  y'acHiOH  and  ./'o'o!  are  multiplied  by  4. 

The  methanol  concentration  distributions  and  the  methanol 
flux  of  these  cases  are  shown  in  Fig.  6.  In  CASE.l  only  a 
small  fraction  of  methanol  penetrates  the  cathodic  catalyst  layer, 
reaches  the  cathodic  flow  channel,  and  escapes  into  the  air.  In 
CASE_0.25,  due  to  the  low  catalyst  efficiency,  a  large  fraction 
of  methanol  escapes.  While  in  CASE_4,  all  methanol  oxidizes 
in  the  cathodic  catalyst  layer.  This  result  indicates  that,  when  the 
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Fig.  5.  Model-predicted  methanol  flux  at  different  cell  current  densities  with 
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catalyst  is  not  efficient,  the  cell  output  voltage  that  is  estimated 
by  assuming  that  the  permeating  methanol  gets  fully  depleted 
in  the  cathode,  as  is  done  in  some  other  models,  may  not  be 
reliable. 

4.  Conclusions 

We  have  developed  a  one-dimensional  mathematical  model 
for  the  DMFC  which  incorporates  fluid  dynamics  calculations  of 
the  flow  channels,  methanol  transport,  oxygen  transport,  reac¬ 
tion  kinetics  within  the  anodic  catalyst  layer  and  the  cathodic 
catalyst  layer  and,  most  important  of  all,  the  mixed  poten¬ 
tial  effect.  The  mixed  potential  was  calculated  by  means  of 
estimating  the  internal  current  generated  by  methanol  perme¬ 
ation  into  the  cathodic  catalyst  layer  without  making  the  usual 
assumption  as  found  in  the  literature  that  methanol  transport  is 
completely  depleted  in  the  cathodic  catalyst  layer.  Predicting 
whether  methanol  can  reach  the  cathodic  flow  channel,  vaporize 
and  exit  the  cell  is  also  made  possible. 

We  have  also  validated  our  model  by  an  experiment  with  a 
home-assembled  cell.  Key  parameter  values  of  the  model  for  this 
particular  cell  are  determined  through  calibration  which  involves 
fitting  of  I-V characteristic  curves  to  experimental  data.  With  key 
parameter  values  known,  we  have  estimated  the  mixed  potential 
effect  and  predicted  methanol  concentration  distributions  and 
methanol  flux  inside  the  cell  at  different  cell  current  densities. 

We  have  also  pointed  out  by  a  case  study  that  predicting  cell 
output  voltage  by  assuming  that  methanol  transfer  gets  depleted 
in  the  cathodic  catalyst  layer  is  not  always  reliable,  especially 
when  the  catalyst  is  not  efficient. 
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Appendix  A.  Fluid  dynamics  for  the  flow  channels 


We  only  consider  laminar  flow  in  a  region  of  the  flow  channel 
where  velocity  and  concentration  are  fully  developed.  The  two- 
dimensional  momentum  equations  are 


P 


du  du  \ 

u - 1-  v  — 

dx  dy ) 


dP  (d2u  d2u\ 
dx  ^  \  dx 2  dy2  ) 


P 


dv  dv\ 
dx  dy ) 


dP  f  d2v  d2v\ 
dy  ^  \  dx2  dy2  ) 


(Al) 

(A2) 


where  Fx  and  F y  are,  respectively,  the  body  forces  acting  on  the 
fluid  in  the  v  and  y  directions,  u  and  v  the  velocity  components  in 
the  v  and  y  directions,  respectively,  P  the  pressure,  p  the  density 
and  /x  the  drag  coefficient. 
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Eq.  (Al)  can  be  neglected  by  applying  u-  0.  Eq.  (A2)  can  be 
simplified  to 


dp  _  /a2 v\ 

a y  ^  \dx2  ) 


(A3) 


by  applying  ^  =  0  and  Fy  =  0.  Assume  ^  =  const,  v(x)  is 
found  to  be 

p(x)  =  ax  +  £>x  +  c, 


where  a,  and  c  are  constants.  With  boundary  conditions  n(0)  = 
0  and  v{l)  =  0,  we  have 


6/ 

d/3 


x(/  —  x), 


(A4) 


where  /  is  the  depth  of  the  flow  channel,  d  the  width  of  the  flow 

channel  and  /  =  d  v  dx  the  flow  rate. 

The  two-dimensional  mass  equation  is 


dc  dc  . 

u  — - 1-  n —  1  =  D 


dx 


dy 


a2c  a2c 
3x2  ay 


(A5) 


where  u  and  n  are  the  velocity  components  in  the  x  and  y 
directions,  respectively,  C  the  concentration  and  D  the  diffusion 
coefficient. 

r)2  f~' 

After  applying  u-  0  and  neglecting  by  assuming  that 

variation  of  C  along  y  is  minimal,  the  mass  equation  is  simplified 
to 


ac  a2c 

v —  =  D — 
dy  dx2 


(A6) 


The  dependences  of  C  on  x  and  y  can  be  isolated  by  assuming 
C(x,y)  =  Cx(x)Cy(y)  and  introducing  a  constant  w  which  is  to  be 
determined  by  boundary  conditions.  Cy(y)  is  found  to  be 


Cy(y)  =  exp  (-y/w), 


(A7) 


and  Eq.  (A6)  becomes 


dfCx 

dx2 


v 

wD 


6/ 

wDdP 


x(l  -  x)Cx. 


(A8) 


Define  Nx  =  —D^,  Eq.  (A8)  becomes 

dNx  6/ 

=  — J-yX(l~x)Cx. 
dx  wdF 


(A9) 


Applying  Eq.  (A9)  to  our  one-dimensional  model  we  get  Eqs. 
(2),  (28)  and  (30). 
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